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In these lectures we introduce the Feynman-Schwinger representation method for solv- 
ing nonperturbative problems in field theory. As an introduction we first give a brief 
overview of integral equations and path integral methods for solving nonperturbative 
problems. Then we discuss the Feynman-Schwinger (FSR) representation method with 
applications to scalar interactions. The FSR approach is a continuum path integral in- 
tegral approach in terms of covariant trajectories of particles. Using the exact results 
provided by the FSR approach we test the reliability of commonly used approximations 
for nonperturbative summation of interactions for few body systems. 

1 Introduction 

Physics research in general is driven by the goal of finding the correct Lagrangian 
for a given system. Once we have a candidate for the correct Lagrangian one must 
be able to relate it to observables. It is clear that the question of "what is the 
correct Lagrangian?" is inseparable from the question of "what is the exact pre- 
diction dictated by a given Lagrangian?" , since we can never be sure that we have 
the right dynamics if we do not know how to calculate the exact result with it. 
However, making exact predictions using a Lagrangian is not always an easy task. 
Therefore, especially in field theory, one has to make approximations. One common 
approximation is known as perturbation theory. Perturbation theory involves mak- 
ing an expansion in the coupling strength of the interaction. It is expected to work 
particularly for small couplings. However, irrespective of how small the coupling 
strength is, it is well known that perturbation theory can not explain bound states. 
This fact can be observed even at the level of classical mechanics. Consider the 
example of a simple harmonic oscillator: In this simple example the Lagrangian is 
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Fig. 1. Simple harmonic oscillator 
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given by 

L = imi 2 - ^gx 2 , (1) 
and from the Euler-Lagrange equations the nonperturbative result follows as: 

x(t) = ASinfwt), w = J — . (2) 

V m 

One might express this result as a power series in the coupling strength 
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(3) 



As this expansion shows a perturbative truncation of the above series can not 
produce a bound state. In other words, in order to be able to obtain a bound state 
result one must sum the interactions to all orders. 

The situation in field theory is similar. Bound states in field theory are identified 
by the pole of Green's function. In general the Green's function can be expanded in 



G = 




Fig. 2. Bound state mass is determined by the pole of Green's function 

powers of the coupling strength. The question is then: Can we make a truncation 
in the perturbation series at order g 2n and still obtain a bound state ? 



G ~ Go + g 2 G 1 +. 9 4 G 2 + • • • + g M G l +■■ ■ + g 2n G n 



(4) 



Since this is a finite series, the bound state singularity could only come from indi- 
vidual terms 
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However this possibility leads to a contradiction since it implies that the dynamics 
of the bound state is independent of g, since g comes only as an overall factor 
in front of an individual term Gi. The case where more than one term is singular 
leads to the same contradiction. Therefore individual G,'s must be nonsingular, and 
bound state singularity must come from an infinite summation of the perturbation 
series. Situation is similar to the expansion 
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where the singularity of the left hand side can not be obtained by a finite truncation 
of the right side. Therefore bound states are always fully nonperturbative. 

With the discovery of quantum chromo-dynamics (QCD) nonperturbative cal- 
culations in field theory have become even more essential. It is known that the 
building blocks of matter, quarks and gluons, only exist in bound states. Therefore 
any reaction that involves quarks will necessarily involve bound states in the ini- 
tial and/or final states. This implies that even at high momentum transfers, where 
QCD is perturbative, formation of quarks into a bound state necessitates a non- 
perturbative treatment. Therefore it is essential to develop new methods for doing 
nonperturbative calculations in field theory. 

The plan of this lecture is as follows. In the following section a review of non- 
perturbative methods in field theory will be given. Later in Sections 3 and 4 the 
Feynman-Schwinger representation will be introduced through examples. In partic- 
ular the emphasis will be on comparing various nonperturbative results obtained 
by different methods. It will be shown with examples that nonperturbative calcu- 
lations are interesting and exact nonperturbative results could significantly differ 
from those obtained by approximate nonperturbative methods. In the last section 
simple perturbation theory results will be derived using the FSR approach. 

2 Nonperturbative methods 

Nonperturbative calculations can be divided into two general categories. These 
are i) Integral equations, and ii) Path integrals. In the following two subsections 
these approaches will be briefly discussed 

2.1 Integral Equations 

Integral equations have been used for a long time to sum interactions to all orders 
with various approximations. [0, ||, |], ||, |) In general a complete solution of field 
theory to all orders can be provided by an infinite set of integral equations relating 
vertices and propagators of the theory to each other. However solving an infinite set 
of equations is beyond our reach and usually integral equations are truncated by 
various assumptions about the interaction kernels and vertices. The most commonly 
used integral equations are those that deal with 1, and 2-body problems. Here we 
give 2 examples for the 2-body bound state problem. The first example is the Bethe- 
Salpeter equation in the ladder approximation, (l) In the ladder approximation the 
Bethe-Salpeter equation sums ladder diagrams to all orders (Fig. ||). Self energies 
vertex corrections and crossed ladder exchanges are left out in this approximation. 
For the simple case of scalar particles, The vertex function &p(k) determining the 
structure of bound state satisfies an integral equation. For the simple case of scalar 
particles bound state equation can be written as 
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Fig. 3. Bethe-Salpeter equation in the ladder approximation sums only ladder type dia- 
grams 

where G(q — k) is the interaction kernel. This equation could be solved using nu- 
merical methods to find the bound state mass P and the vertex function <3>p(fc). 
The vertex function is similar to the quantum mechanical wave function. In prin- 
ciple it contains all information about the bound state. A serious deficiency of the 
ladder approximation is that it does not have the correct one body limit when 
one of the particles is infinitely heavy. || This is due to the fact that the ladder 
approximation ignores the crossed ladder type exchanges between the particles. As 
it will be shown by explicit examples in this article the crossed ladder exchanges 
play a crucial role in obtaining the correct result not only in the 2-body problem 
but also in the calculation of self energies for the 1-body propagators. 

The problem of 1-body limit when one of the particles is infinitely heavy is solved 
by the Gross equation (Fig. ||). || In the Gross equation the heavier constituent is 
constraint to its physical mass shell (q 2 = M 2 ) Putting the heavier constituent on 





Fig. 4. In the limit of infinitely heavy-light systems Gross equation effectively sums all 
ladder and crossed ladder diagrams. 
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its mass shell and summing only ladder diagrams effectively is equivalent to sum- 
ming all ladder and crossed ladder diagrams. When the heavy particle is constraint 
to its mass shell, the bound state equation takes the following form 



When bound state of equal (or close) mass particles are under consideration the 
Gross equation can be symmetrized by picking up the mass pole contribution of 
both particles in different channels. Gross equation is a manifestly covariant rela- 
tivistic equation, and in the nonrelativistic limit Schroedinger equation is recovered. 
In the literature this equation has been successfully used to analyze relativistic 
bound states. 

Some general comments on integral equations are as follows: There are a number 
of useful features of integral equations. Integral equations are a natural extension 
of nonrelativistic quantum mechanics to field theory. They are a practical tool for 
modeling and doing simple calculations using the tools of field theory. Because of 
their similarity to Schroedinger equation in quantum mechanics, integral equations 
in field theory provide an intuitively clear picture of physics. Furthermore the nu- 
merical cost of solving integral equations is negligibly small compared to methods 
such as path integrals (to be discussed in the next section). The last important 
advantage of integral equations is the fact that they can be solved in Minkowski 
metric. The only other alternative to integral equations, path integrals, make use 
of the Euclidean metric which limits their applicability. This is an important tech- 
nical problem -particularly for calculations of scattering reactions, form factors and 
decays- since physical particles live in the timelike region. 

On the other hand integral equations have some drawbacks. The first problem is 
that they are not exact. Without knowing the exact result it is not possible to claim 
that integral equations are even a good approximation to the full theory. The second 
problem is that integral equations in general do not respect the symmetries of the 
underlying Lagrangian. Except for very special approximations, integral equations 
are not gauge invariant. Therefore a more rigorous and systematic approach is 
needed. This is where path integrals play a significant role. In the next section we 
introduce the method of path integrals in quantum mechanics and field theory. 

2.2 Path Integrals 

Path integrals provide a systematic method for summing interactions to all or- 
ders. First we begin by giving path integral expressions for quantum mechanics. 
The matrix element for a transition from an initial state \qi,U) to a final state 
!?/>*/) is S iven b y 





(qf,tf\qi,ti) H 



\<4ti H/Si 



(9) 




(10) 
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where subscripts H, and S refer to Heisenberg and Schroedinger states. In this 
path integral expression all trajectories contribute to the final result with equal 
weight. Classical limit can be readily obtained by letting fi — > 0. In this limit small 




Fig. 5. In quantum mechanics all paths contribute. 



variations in the action will lead to large oscillations and the dominant contribution 
will come from those trajectories that minimize the action: 



SS 



■is 



L[q,q,t]dt = 0. 



(11) 



The minimization of the action leads to Euler-Lagrange equations of classical me- 
chanics: 



d 
dt 



dL 
dq 



dL 



0. 



(12) 



Keeping the field theory applications in mind we consider the time ordered products 
of operators. The time ordered product of operators can also be expressed in the 
form of a path integral: 



(qf, t f\ T (Q( t i)q(t 2 )---q(t n )\qt,t l ) = \Vqq(t 1 )q(t 2 )---q{t n )exp 



dtL{q, q) 



In field theory applications the ground state expectation values of time ordered 
products are particularly needed. In order to show how the ground state expectation 
value can be obtain let us introduce a complete set of energy eigenstates into the 
2-point Green's function 



(g/,t|T[$(ti)$(t 2 )]|<?i,-t) = Yl <?/.*/|n><n| T [5(*i)$(t2)]|m)<m|ft,ti).(13) 
Noting that 



n,m— 



(q,t\n) 
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(14) 
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one obtains 

oo 

{q f ,t\TW±W2)]\qu-t) = J2 c- l(£ " +Bm) V™(9/)V'mfe)^|T[ ( z(ti)g(i2)]|m). 

n,m— 

As t — > — ioo the ground state (n = 0) will give the dominant contribution, 



lim (g/,f|T[g(*i)g(t 2 )]| gi ,-t) 



-2_E t 



V'o(<7/)V'o(ft)(0|T[g(t 1 )g(t 2 )]|0).(15) 



<<//,*ki,-*> 

Therefore the ground state expectation value is obtained as 

(qf,°°\T[q(ti)q(t 2 )}\qi,-oo) 



(q f ,oc\qi,-oo) 



(16) 



|©««(tl)?(t2)"- 


g(i„) exp 


J — oo 




J exp 


i J x C o dtL ( ( i><i) 





In general the ground state expectation value of the time ordered products of an 
arbitrary number of operators can be written as 



(0\T[q(t 1 )q(t 2 )---q(t n )}\0) 



In going from quantum mechanics to quantum field theory we recognize that a field 
is an infinite array of coordinates, such as an infinite array of coupled oscillators, 
and particles are associated with the normal modes of oscillators. Therefore the 
quantum mechanical derivation can be generalized by replacing 

q{t) -►^(xi,t),^(x2,t),"-,0(x„,t),-" ) (17) 

where <j>(xi,t) is the displacement of the oscillator at Xi. The generalized path 



<KM) (KM 



♦CM) 



Fig. 6. Field as the coordinates of an infinite array of coupled oscillators 

integral in terms of fields is given by: 

-iH(t f -ti) |0(x,ti)) = jp^exp 

A 



(<j)(x, i/)|exp 
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where 



(19) 



represents a sum over all possible field configurations. The step of going from parti- 
cle trajectories q(t) to fields 0(x, t) dramatically increases the dimensionality of the 
problem. While quantum mechanical path integral sums all possible trajectories 
(line configuration) , field theoretical path integral sums all possible field configura- 
tions (volume configuration). 

The Green's function in field theory is given by the path integral expression: 



(0\T[<t>( Xl )<l>(x2) ■ ■ ■ <l>(x T 



This result can be obtained from a generating function 

6 n Z[J(x)] 



j[2ty]#sl)0(ta)-" 


4>{x n ) exp 


i J d x L(x) 


j [D(j>] exp 


i J d 4 x L(x) 





(O\T(<j)(x 1 )^x 2 )---0(x n )\O) 



5J(X\)8J(X2) ■ ■ ■ 5J{x n ) 



J=0 



where 



Z[J(x)] = [D$exp 



| [V<j>] exp i J 



i / d 4 x [L(x) + J(x)<f»(x)] 



(20) 



(21) 



While path integrals provide a compact expression for the exact nonperturbative 
result for propagators, evaluation of the path integral is a nontrivial task. In general, 
field theoretical path integrals must be evaluated by numerical integration methods, 
such as Monte-Carlo integration. The best known numerical integration method is 
lattice gauge theory. || Lattice gauge theory involves a discretization of space-time. 
In a discretized lattice particles can be located on the discrete lattice sites and ex- 











r 





Fig. 7. In lattice gauge theory space-time is discretized within a finite box. 
change fields are represented by links between the sites. After discretization the 
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path integral can be performed essentially by a brute force method. Particularly 
for QCD lattice gauge theory is currently the only method that can produce non- 
perturbative results starting directly from the QCD Lagrangian. However lattice 
calculations are not without drawbacks. Discretization of space-time by a cubic 
lattice violates rotational symmetry. In addition the cost of computations critically 
depend on the size of the lattice. Because of this limitation complex applications 
such as calculation of form factors or scattering reactions are beyond the reach of 
current lattice applications. Finally matter loops are not accounted for (quenched 
approximation) in most lattice calculations due to high cost. 

In the next section we present a more efficient method of performing path in- 
tegrals in field theory for simple scalar interactions. This method is known as 
Feynman-Schwinger representation(FSR). g g §, 0, |H H H @ Through 
applications of the FSR, the importantance of exact nonperturbative calculations 
will be shown with explicit examples. 



3 Feynman-Schwinger representation approach 

The basic idea in the FSR approach is to transform the field theoretical path 
integral such that a quantum mechanical path integral in terms of particle trajec- 
tories is obtained. In this section we consider application of the FSR technique to 
scalar QED. The Minkowski metric expression for the scalar QED Lagrangian in 
Stueckelberg form is given by 



C 



SQED 



-m 2 Y 2 - -F 2 
4 



1 



2 A 2 



2 F 



where A represents the gauge field of mass and x is the charged field of mass 
m. The presence of a mass term for the exchange field breaks the gauge invariance. 
Here the mass term was introduced in order to avoid infrared singularities when 
application to 0+1 dimension is considered later in the next section. For dimensions 
larger than n=2 the infrared singularity does not exist and therefore the limit /i — > 
can be safely taken to insure gauge invariance. 

The path integral is to be performed in Euclidean metric. Therefore we perform 
a Wick rotation: 



exp 



d 4 xC 



M 



exp 



J d i x Ce 



The Wick rotation for coordinates is obtained by 

-ix 
id . 



5n = 



x 
d 



dx 



(22) 
(23) 



The transformation of field A under Wick rotation is found by noting that under a 
gauge transformation: 



A, 



A, 



(24) 
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Then, under a Wick rotation: 

A iA . 
The Wick rotated Lagrangian for SQED is given by: 

Csqed = X m 2 -8 2 -2ieA8-iedA + e 2 A 2 X + £a- 
The exchange field part of the Lagrangian is given by 

C A = \A^ 2 g^-\d^d v )A v + - A F 2 , 

= ^[(m 2 - U) 9llv + (1 - X)d IM d v }A v . 
We employ the Feynman gauge A = 1 which yields 

C A = \A v {p 2 - U)A V . 



(25) 
(26) 

(27) 
(28) 

(29) 



The two-body Green's function for the transition from an initial state $i to final 
state $/ is given by 



G(y,y\x,x)=N J V X * J V X J VA 



-Se 



where 



Se — J d A x Csqed, 
and a gauge invariant 2-body state <I> is defined by 

$(x,x) = x*(x)U(x,x)x(x). 



(30) 
(31) 

(32) 



The gauge link U(x,y) which insures gauge invariance of bilinear product of fields 
is defined by 



— ie / dz A{z) 



U (x, y) = exp 



One can easily see that under a local gauge transformation 



eA(x] 



X(x) 



<&i(x,x) remains gauge invariant 



$(x, x) — > exp 



-ieA(x) + ieA(x) — ie dz^ d^A 



r 



$(x,x). 
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(34) 
(35) 



X*(x)U(x,x)x(x) (36) 
(37) 
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Performing path integrals over \ an d X* fields in Eq. [50] one finds 

G(y,y\x,x) = N^VA (detS) U(x, x)U*(y, y)S(x, y)S(x, y) e" 5 ^, (38) 
where the interacting 1-body propagator S(x, y) is defined by 

s ( x >v) = (y\ 2-Afp * I*). ( 39 ) 

+ H(z,p) 

H(z,p) ee {p + ieA{z)f. (40) 

The Green's function Eq. ^8| in principle includes contributions coming from all 
possible interactions. The determinant in Eq. ^8| accounts for matter loops and in 
the quenched approximation it is set equal to one (detS* — > 1). Analytical calculation 




Fig. 8. Various interactions in presence of a matter loop 



of the path integral over gauge field A in Eq. |3§| seems difficult due to nontrivial 
A dependence in S(x,y). In more complicated theories such as QCD integration of 
gauge field integral, as far as we know, is not analytically doable. Therefore, in QCD, 
the only option is to do the gauge field path integral by using a brute force method 
on a discretized space-time lattice. However for the simple scalar QED interaction 
under consideration it is in fact possible to go further and eliminate the path integral 
over field A. In order to be able to carry out the remaining path integral over 
the exchange field A it is desirable to represent the interacting propagator in the 
form of an exponential. This can be achieved by using a Feynman representation 
for the interacting propagator. The first step involves the exponentiation of the 
denominator in Eq. g9f 

S(x,y) = / ds e- sm2 (y|exp[-,stf]|x), (41) 
Jo 

This expression is similar to a quantum mechanical propagator where s = it, and 
H is a covariant Hamiltonian in terms of 4-vector momentum and coordinates. It 
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is known that one can use a path integral representation for quantum mechanical 
propagators. A covariant Lagrangian can easily be obtained from the Hamiltonian 

z 2 

H(z,p) = (p + ieA(z)r, L(z,z) = — - iezA(z). (42) 

Using this Lagrangian a path integral representation for the interacting propagator 
can be constructed 



I>1 



S(x, y) = ! ds Vz exp 



2 1 

-sm 



dr z 2 {t) — ie I dr zA(z(t)) 



(43) 



where the boundary conditions are given by z(0) = x, z(s) = y. This representa- 
tion allows one to perform the remaining path integral over the exchange field A. 
The final result for the two-body propagator involves a quantum mechanical path 
integral that sums up contributions coming from all possible trajectories of particles 

G = | ds | ds j(Vz) xy j(Vz) S ye- K ^- K ^{W(C)), (44) 

where the kinetic term is defined by 

1 f 1 

K[z,s] = m 2 s + — / dri 2 (r), 
4s Jo 

and the Wilson loop average (W(C)) is given by 

(W(C)) = J VAexp -ie <j> dz A(z) -^J d 4 z A{z){^ 2 - d 2 )A{z) 



(45) 



(46) 



where the contour of integration C (Fig. ^) follows a clockwise trajectory x — * y — > 
y — > x — > a; as parameters r, and f are varied from to 1. The A integration in 




Fig. 9. Wilson loop. 
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the definition of the Wilson loop average is of standard gaussian form and can be 
easily performed to obtain 



(W(C)) = exp 



/ dzp I dz v A /J , JJ (z - z, (i) 
Jc Jc 

e ipx 



( dp e 
9]XV ) J2^p^~ 



(47) 
(48) 



The self energy and the exchange interaction contributions, which are embedded in 
expression |47], have different signs. This follows from the fact that particles forming 
the two body bound state carry opposite charges. By the result given in Eq. [44] 
and Eq. ^ path integration expression involving fields has been transformed into 
a path integral representation involving trajectories of particles. The bound state 
spectrum can be determined from the spectral decomposition of the two body 
Green's function 

oo 

G(T) = 5> n e-" l " T , (49) 

71 = 

where T is defined as the average time between the initial and final states 

T = 7}(V4 + Vi ~ x 4 - x 4 ). (50) 

In the limit of large T, the ground state mass is given by 

fT>ZS'[Z]e- s W 



l -^oo al 



3.1 Application to 0+1 dimensions 



JVZe- s W 



(51) 



Massive scalar QED in 0+1 dimension is a simple interaction that enables one 
to obtain a fully analytical result for the dressed and bound state masses within 
the FSR approach. In this section we compare the self energy result obtained by 
approximate methods with the full result obtained from the Feynman-Schwinger 
representation. In space + 1 time dimension there is no continuum spectrum, and 
only bound states exist. 

In general Wilson loop average depends on the trajectory of particles: 



(W{C)) 



exp 



dz dzA(z — z,fi) 
Jc Jc 



(52) 



In 0+1 dimensions Wilson loop integral is essentially a line integral and all trajec- 
tories contribute equally. Therefore in 0+1 dimension Wilson loop average does not 
depend on the shape of trajectory 



(W{C)) = exp 



f T f T 

— / dz dzA(z — z,/i) 

2 Jo Jo 



(53) 
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Notice that in 0+ld all trajectories lie on a straight line. Possible variations can only 
come from those trajectories which fold onto themselves. However the contribution 
of the folded sections of trajectories identically vanish. In addition, contribution 
of matter loops to the Wilson loop average is identically zero. The typical loop 
contribution in 1-dimension can be written as 



T 



dz( dz' + dz')A(z-z') = 0. (54) 



The vanishing of matter loop contribution implies that the quenched calculations 
give exact results. Therefore massive SQED in 0+1 dimension is a remarkably sim- 
ple example where we can compare exact analytic solutions of field theory with 
various approximate nonperturbative methods. Furthermore, analytical results in 
0+ld provide a test case for the numerical routines which are used in higher di- 
mensions. 

In 0+ld the interaction kernel A(z — z, fi) is given by 

e -/j,\z-z\ 

A(z-z,M) = 2 • (55) 

Using this kernel the Wilson loop average Eq. [53| is calculated exactly: 

e 2 T / 1 - e-» 



(W(C)) = exp 



1 - 



2/i 2 V ptT 



(56) 



Remaining integrals over ds and Vz in Eq. |4| provide the free particle exponential 
fall of e~ mT at large times. Therefore the spectrum of the Green's function can be 
trivially calculated using Eqs. |4], [5l], and |5(| Exact analytic FSR results for bound 
state masses are as follows: 

e 2 

1-body: Mi=m+— , (57) 

2 - body : M 2 = 2m, (58) 

e 2 

n — body : M n = nm — n (n — 2) — . (59) 

Bound state mass results given above include the self energy contributions. The 
n-body result shows that for any given coupling strength e, there is an upper limit 
in particle number n beyond which results become unstable. Dynamics of the two 
body bound state mass generation is similar to pion mass generation. Self energies 
of particles exactly cancel the binding energy to produce a bound state mass that 
is proportional to the current masses 

2 2 2 

(to + — ) + (to + — ) = 2m. (60) 

Exact massive SQED results share some common features with QCD in 1+1 d 113] , 
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— /i — > ==>■ Mi — > oo, where plays the role of an infrared cut-off. 

— 2-body bound state mass M2 is independent of ji, 

— When "cut-off" /1 — > =>■ no bound states for n > 2. 

— When current masses vanishes m — > ==> M2 = "chiral limit" 

This toy model provides a possible test case also for the lattice gauge theory cal- 
culations. 

3.2 Comparison of the exact FSR results with approximate nonperturbative 
methods 

In this section we take a closer look at 1-body mass pole calculations. Popular 
methods frequently used in finding the dressed mass of a particle is to do a sim- 
ple bubble summation or solve the 1-body Dyson-Schwinger equation in rainbow 
approximation. It is interesting to compare results given by the bubble summation 
and the Dyson-Schwinger with the exact FSR result. Below we first give a quick 
overview of how dressed masses can be obtained in bubble summation and the 
Dyson-Schwinger equation approaches. 



Bubble sum: 



• * • 



Dyson-Schwinger: 



Feynman-Schwinger: 



Fig. 10. Various interactions included in each approach are shown. 

The simple bubble summation involves a summation of all bubble diagrams 
(Fig. [n]) to all orders. The dressed propagator is given by 

Mp) = 2 ^ 2.WV (61) 
The dressed mass M is determined from the self energy using 



M = yjm 2 + E(iAf). (62) 
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The self energy for the simple bubble sum is given by 



2 f°° dk 1 / (2p-k) 2 



The self energy integral in this case is trivial and can be performed analytically, 
and the dressed mass is determined from Eq. |6^ 

The rainbow Dyson-Schwinger equation sums more diagrams than the simple 
bubble summation (Fig. |To| ) . The self energy of the rainbow Dyson-Schwinger equa- 
tion involves a momentum dependent mass. 



dk 1 f (2p-k) 2 
, fcr (k 2 + fj, 2 ) \ [(p - fc)2 + m 2 + E(jfc)] 



S(P) - ^ „, 2 , „ 2 J FT- ,, , . ,„2 . ^TTTT - 1 • (64) 



In this case the self energy is nontrivial and it must be determined by a numerical 
solution of Eq. ^54|. The dressed mass is determined by the logarithmic derivative of 
the dressed propagator in coordinate space 

M = - lim Al 0g [A d (*)]. (65) 
T— >oo al 




e 2 (GeV 3 ) 



Fig. 11. The function M(e 2 ) calculated by the FSR approach, the Dyson-Schwinger 
equation in the rainbow approximation, and the bubble summation for values of m = n = 1 
GeV. While the exact result is always real, the rainbow DSE and the bubble summation 
results become complex beyond a critical coupling. 
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The type of diagrams summed by each method is shown Fig. [10[ Note that the 
matter loops do not give any contribution as explained earlier. Results obtained 
by these three methods are shown in Fig. [ll]. It is interesting to note that the 
simple bubble summation and the rainbow Dyson-Schwinger results display similar 
behavior. While the exact result provided by the FSR linearly increases for all 
coupling strengths, both the simple bubble summation and the rainbow Dyson- 
Schwinger results come to a critical point beyond which solutions for the dressed 
masses become complex. This example very clearly shows that conclusions about the 
mass poles of propagators based on approximate methods such as rainbow Dyson- 
Schwinger equation can be misleading. 
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Fig. 12. As a test of the numerical methods the numerical Monte-Carlo result for 1-body 
dressed mass in SQED is compared with the analytical result. 

The simplicity of massive SQED in 0+ld also provides an excellent opportunity 
to test the numerical Monte-Carlo integration methods which are normally needed 
at higher dimensions. Therefore as a test of the numerical methods used for the 
remainder of applications in this work, we demonstrate in Fig. 12 that the numerical 
Monte-Carlo results for 1-body dressed mass in SQED correctly reproduce the 
analytical result. In Fig. [l^ the time dependence of the dressed mass given by 
Eq. [5l] is shown. Parameters used for this plot are /1 = 0.15 Gev, and e = 0.15 
GeV 15 . 

In the next section we consider the application of the FSR approach to scalar 
X </> interaction 
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4 Scalar x 2 4> interaction with the FSR approach 



We consider the theory of charged scalar particles X °f mass m interacting 
through the exchange of a neutral scalar particle <f> of mass fi. The Euclidean La- 
grangian for this theory is given by 



(66) 



The 2-body propagator for the transition from the initial state $i = x*{ x )xi x ) to 
final state $/ = X*{v)x{v) is given by 



d 4 xC 



(67) 



G(y,y\x,x)=N\Vx*\V X \VA exp 

After the usual integration of matter fields is done the Green's function reduces to 

G(y,y\x,x) = N^Vcj> (det£) S(x, y)S(x, y) e~ 8 ^ . (68) 

As in the case of scalar QED we employ the quenched approximation: detS — > 1. 
The interacting propagator S(x, y) is defined as 

1 



S(x,y) 
H(z,p) 



m? + H(z,p) 



\x) 



P 2 -gff>{z). 



(69) 
(70) 



We exponentiate the denominator by introducing an s integration along the imag- 
inary axis with an e prescription 



S(x,y) = 



ds e - s(m2+ie) (t/|exp[-sff] \x). 



(71) 



This representation should be compared with the representation used earlier in 
SQED Eq. |d[ Here the integration is done along the imaginary axis because H is 
not positive definite. Again a quantum mechanical path integral representation can 
be constructed by recognizing that Lagrangian corresponding to H Eq. is given 
by 



L{z,z) = — + g<j){z). 
The path integral representation for the interacting propagator is 



(72) 



S{x,y) 



'}>} 



T>z exp 



is(m 2 + ie) — — / dTZ 2 (r)+ig I dr 4>(z(t)) 



This representation allows the elimination of the integral over exchange field <j>. The 
2-body propagator reduces to 



G = - 
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where mass and kinetic term is given by 

1 r 1 

K[z,s] = (m 2 + ie)s-—j dr z 2 (t). 
The field integration 1$ is a standard gaussian integration 

I* = jp^exp +ig(J dr0(z(r))+ J dfcj>{z(f))} S'n] 
= expf -Vo[z, s] - 2Vi2[z,z,s,s] - V [z,s] I, 



(74) 

(75) 
(76) 



where Vq and V\2 (self and exchange energy contributions in Fig. |T^) are defined 
by 



V [z,s] 



2 



Vi 2 [z,z,s,s] = — ss 



{ dr { dr' 
J o Jo 

Jo Jo 



A(z(t)-z(t'),m), 



df A(z(r)-f(f),M). 



(77) 
(78) 



It should be noted that the interaction terms explicitly depend on the s variable, 
which was not the case for SQED. The interaction kernel A is given by 



d 4 p e ip - a 



(79) 



(2tt) 4 p 2 + M 2 47r 2 |a;|' 
In order to be able to compute the path integral over trajectories a discretization 



z(T) 




Fig. 13. Sample trajectories with self and exchange interactions, 
of the path integral is needed 

(Vz) xy (N/4-ks) 2N U^^ 
Czech. J. Phys. 50 (2000) A 



d A Zi 



(80) 
19 



Qetin §avkh 



The s dependence is crucial for correct normalization. After discretization the 1- 
body propagator takes the following form 



G 



im s 



4s 



s 2 v 



(81) 



This is an oscillatory and regular integral and it is not convenient for Monte-Carlo 
integration. The origin of the oscillation is the fact that s integral was introduced 
along the imaginary axis, 



Rep. 1 : S(x,y) =<y\ 



1, 



ds exp 



— s(m 2 — d 2 + g4> + ie) 



| x > . (82) 



In earlier works [TlJ a nonoscillatory Feynman-Schwinger representation was 
used, 



Rep. 2 : S(x, y) =< y\\ ds exp 



-s(ra 2 - d 2 + g4>) 



x > 



Rep. 2 leads to a nonoscillatory and divergent result 

ds 2 k 2 

Z2N eX P - m S - 

o " L 



G 



i; 



4s 



+ s 2 v 



(83) 



(84) 



and the large s divergence was regulated by a cut-off A. This is not a satisfactory 
prescription since it relies on an arbitrary cut-off. Later it was shown [[l3[ [l5| that 
the correct procedure is to start with Rep.l and make a Wick rotation such that 
the final result is nonoscillatory and regular. The implementation of Wick rotation 
however is nontrivial. Consider the s-dependent part of the integral for the 1-body 
propagator 



G oc 



i; 



ds 

72lv ex P 



k 2 

2 - 2 

im s — i s v 

4s 



(85) 



It is clear that a replacement of s — > is leads to a divergent result. The problem with 
Wick rotation (Fig [ti]) comes from the fact that the s integral is infinite both along 
the imaginary axis and along the contour at infinity. These two infinities cancel to 
yield a finite integral along the real axis. As g — > the dominant contribution to 
the s integral in Eq. comes from the stationary point 



is 



■ k 

2m 



(86) 



Therefore one might suppress the integrand away from the stationary point by 
introducing a damping factor R 



20 



R(s,s ) 
A 



i - (s - is ) 2 /r 2 



(87) 
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Im(s) 




Fig. 14. Wick rotation in s integration 



With this factor the integral in Eq. 

/ oo , 

I as 



G 



a is modified as 
k 2 



C 2N 



exp 



im 2 s 



4s R 2 (s,s )v 



(88) 



This modification allows us to make a Wick rotation since the contribution of the 
contour at infinity now vanishes. However this procedure relies on the fact that 
there exist a stationary point. It can be seen from the original expression Eq. |85| 
that this is not always true. According to the original integral the stationary point 
is given by the following equation 



im +i- 



4s 2 



2sv = 0. 



(89) 



The stationary point is determined by the first intersection of a cubic plot with 
the positive s axis as shown in Fig. [l^. The plot in Fig. [l5] shows that as coupling 
strength is increased the curve no longer crosses the positive s axis. Therefore be- 
yond a critical coupling strength the stationary point vanishes and mass results 
should be unstable. Having noted that the original expression Eq. |85| has a crit- 
ical point, we now turn to perform a Wick rotation on the modified expression 
Eq. ^8|. Wick rotation in Eq. ^8] amounts to a simple replacement s — > is, and a 
nonoscillatory and regular integral is found: 



G 



i; 



ds 

75jv ex P 



-m 2 s — 



4s R 2 (is,s ) 



(90) 



At first look it seems that the new integral always has a stationary point determined 
by the following equation 



— m + 



4s 2 



2sv 



1 



i? 2 (is,s ) 



O 

S V 



(R 2 (™,s ))' 

R 2 (is,s ) 



= 0. 



(91) 



The key point to remember is that the stationary point we find after the Wick 
rotation should be the same stationary point we had before the Wick rotation. 
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Fig. 15. As coupling strength is increased the stationary point disappears. 

This is required to make sure that the physics remains the same after the Wick 
rotation. Therefore self consistency requires that the stationary point after the Wick 
rotation is at s — isg. In that case R(iso, sq) = 1, and (i? 2 (iso, so))' = and the 
equation determining the critical point Eq. ^l] reduces to the earlier original form 
given by Eq. |8^. Therefore self consistency requirement guarantees that the critical 
point still exists after the Wick rotation. 

The regularization of the ultraviolet singularities are done using Pauli-Villars 
regularization prescription. Pauli-Villars regularization is particularly convenient 
for numerical integration since it only involves a change in the interaction kernel 

A(x,fi) — ► A(x,n)- A(x,an). (92) 

4.1 Numerical applications of y 2 ^ interaction 

Applications of x 2 </> interaction in 3+ld require numerical Monte-Carlo inte- 
gration. First step is the discretization of particle trajectories, where boundary 




Fig. 16. Number of steps a particle takes between initial and final coordinates is dis- 
cretized. The space-time is continuous and there are no space-time boundaries 

conditions are given by 

2o = x = (xi,x 2) x 3 ,0) z N = y = (2/1,2/2, J/3. T). (93) 
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Discretization employed in the FSR is for the number of steps a particle takes 
between the initial and final states in a 4-d coordinate space. This is very different 
from the discretization employed in lattice gauge theory. Contrary to lattice gauge 
theory, in the FSR approach space-time is continuous and the rotational symmetry 
is respected. An additional benefit is the lack of space-time boundary. This is an 
important advantage of the FSR approach. The lack of space-time lattice boundary 
allows analysis of arbitrarily large systems using the FSR approach. This feature 
provides an opportunity for doing complex applications such as calculation of form 
factors using the FSR approach. 
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Fig. 17. The dependence of the peak of the s-distribution on the coupling strength is 
shown. The peak location is given by s — CT /2m. Beyond the critical coupling strength 
of g 2 — 31GeV 2 a self consistent determination of C is not possible. Therefore beyond the 
critical coupling strength 1-body mass becomes unstable. 

In doing Monte-Carlo sampling we sample trajectories (lines) rather than gauge 
field configurations (in a volume). This leads to a significant reduction in the nu- 
merical cost. The ground state mass of the Green's function is obtained using 

fVZS'[Z]e- s W 

m ° = fvze-sm ■ (94) 

Sampling of trajectories is done using the standard Metropolis algorithm. Metropo- 
lis algorithm insures that configurations sampled are distributed according to the 
weight e~ s I z l. In sampling trajectories the final state (spacial) coordinates of par- 
ticles are integrated out. Integration over final state coordinates puts the system 
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at rest and projects out the s-wave ground state. As trajectories of particles are 
sampled wave function of the system can be determined simply by storing the final 
state configurations of particles in a histogram. 
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Fig. 18. The coupling constant dependence of the 1-body dressed mass is shown. Beyond 
the critical coupling strength of g 2 — 31 GeV 2 the 1-body mass becomes unstable. The 
perturbative bubble summation also has a critical point near g 2 = 490 GeV 2 . 

In sampling trajectories the first step is termalization. In order to insure that 
the initial configuration of trajectories has no effect on results initial updates are 
not taken into account. Depending on the dimensionality of the problem and the 
coupling strength the number of initial updates neeeded for termalization is of the 
order of 1000 updates. In order to satisfy self consistency regarding the location 
of the stationary point discussed earlier, the location of the stationary point must be 
determined carefully. The stationary point can be parametrized as sq = CT/(2m), 
where T/2m is the location of the stationary point when the coupling strength g 
goes to zero. As the coupling strength is increased the stationary point moves out 
(see Figs. [f5[ and |l7| ) and eventually the critical point is reached beyond which 
there is no stationary point. In order to be able to do Monte-Carlo integrations an 
initial guess must be made for the location of the stationary point. Self consistency is 
realized by insuring that the peak location of the s distribution in the Monte-Carlo 
integration agrees with the initial guess for the stationary point. |l5| In Fig. |lj] 
the dependence of the location of the stationary point on the coupling strength is 
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shown. Fig. |l7| shows that beyond the critical point g 2 ~ 31 GeV 2 C goes to infinity 
implying that there is no stationary point. A similar critical behavior was also 
observed in Refs. |^(| within the context of a variational approach. In Fig. [l8| exact 
1-body dressed mass results are shown for m x — 1 GeV, [1$ = 0.15 GeV. Results 
indicate that the perturbative bubble summation deviates from the quenched FSR 
result very significantly. These results are all for a Pauli-Villars mass of 3/z. In 
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Fig. 19. The coupling constant dependence of the 2-body bound state mass is shown. 
Beyond the critical coupling strength of g 2 = 100 GeV 2 the 2-body mass becomes unstable. 
The Bethe-Salpeter equation in the ladder approximation gives the lowest binding. 

Figure |l9| we present the comparison of the 2-body bound state masses obtained 
by the FSR with various bound state equations. The FSR calculation involves 
summation of all ladder and crossed ladder diagrams, and excludes the self energy 
contributions. According to Figure p~s| all bound state equations underbind. Among 
the manifestly covariant equations the Gross equation gives the closest result to the 
exact calculation obtained by the FSR method. This is due to the fact that in the 
limit of infinitely heavy-light systems the Gross equation effectively sums all ladder 
and crossed ladder diagrams. Equal-time equation also produces a strong binding 
but the inclusion of retardation effects pushes the Equal-time results away from the 
exact results (Mandelzweig- Wallace equation JlG|]). In particular the Bethe-Salpeter 
equation in the ladder approximation (BSE in Figure p"9| ) gives the lowest binding. 
Similarly the Blankenbecler-Sugar-Logunov-Tavkhelidze equation JI7L [l8| (BSLT) 
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gives a very low binding. A comparison of the ladder Bethe-Salpeter, Gross, and 
the FSR results shows that the exchange of crossed ladder diagrams plays a crucial 
role. 
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Fig. 20. 3-body bound state results for 3 equal mass particles of mass 1 GeV. 

In Figure |2(] the 3-body bound state results for 3 equal mass particles of mass 
1 GeV is shown. For 3-body case the only available results are the Schroedinger 
and Gross equation results. According to results presented in Figure ^ bound state 
equations underbind for the 3-body case too. Gross equation gives the closest result 
to the exact FSR result. Determination of the wavefunction of bound states is done 
by keeping the final state configurations of particles in a histogram. For example, 
for a 3-body bound state system, the probability distribution of the third particle 
for a given configuration of first and second particles is shown in Fig. In the first 
plot of Fig. |l] two fixed particles are very close to each other such that the third 
particle sees them as a point particle. However as the fixed particles are separated 
from each other the third particle starts having a nonzero probability of being in 
between the two fixed particles (second and third plots of Fig. ^l|). Eventually 
when the two fixed particles are kept away from each other the third particle has a 
nonzero probability distribution only at the origin (the last plot shown in Fig. |^) . 
Until this point the FSR method has been derived and various applications to 
nonperturbative problems have been presented. In the next section, as a check of 
the FSR method, we will obtain perturbative results using the FSR. 
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Fig. 21. Evolution of the probability distribution for the 3rd particle is shown as the 
distance between the two fixed particles is increased. When the fixed particles are very 
close to each other the third particle sees them as a point particle (the upper left plot). 
As the fixed particles are separated from each other the third particle starts penetrating 
between them (2nd and 3rd plots) , and as the two fixed particles are maximally separated 
the third particle spends most of its time in between the two fixed particles (the lower 

right plot). 



5 Perturbative expansion from the FSR approach 

In this section we will show that perturbation theory results can be obtained 
from the FSR expressions. Let us first consider the perturbative expansion of the 
1-body self energy in \ 2 (f> interaction. The exact Greens function can be expanded 
in a power series in g 2 : 

G | ds jpzexp[-m 2 s - drz 2 + y J dr J dr' A(z(t) - z(t')] 

= G (0) +G (1) + •••. (95) 
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The leading contribution (O(g )) is given by: 



G (0) = I dsjvze^[-m 2 s-^j^ drz 2 ], 

(°° ( N \ 2N N N 

= J Hw n^ 4 z„exp[-m 2 S --g 



(zt-Zi-!) 2 ]. (96) 



Noting that: 

^dzidz2 ■ ■ ■ dz n exp 

1 V2 



-A{(«i - x) 2 + (z 2 - zi) 2 + • • • + (z„ - y) 2 } 



(n+ 1)A™ 

all z integrals can be performed and the free (0(5°)) propagator is found as 

G<°> = I ^exphm 2 ,--^-^] (97) 
m 



✓ OC 

J 



Ait 2 \x — y 



■Kiimlx - y\), 



(98) 



which is the free propagator for a massive scalar particle in 3+ Id. Now let us 
consider the next to leading order 0(g 2 ) contribution to the 1-body propagator. 
Order 0(g 2 ) term is given by 

G (1) = dsivzi dr \ dr'exp[-TO 2 s- - f drz 2 ] A(z(r) - z{t')). 

2 J J Jo Jo 4 ^o 



This expression has the following structure 

„2 / 00 



2 /• 00 ✓ s ✓ s 

G« = |-J dsjdTjdT'f(s,T,r'). 



Using the identity 



/ OO / S / OO / OO 

i ds\ drg(s,T)= ds\ drg(s + T,T 
J J Jo Jo 



(99) 



(100) 



one may write 



y- OO S S S T S OO S OO S OO 

j (is dr dr' f{s,T,r') = \ ds\ dr dr' f(s + r + r' \t + t' ,t'). 
J0J0J0 J0J0J0 

Therefore limits of all integrals in can be extended to infinity 

/CO /OO s OO 

G (1) = g 2 \ ds\ drl dT'[f(s + T + T',T + T',T') (101) 
Jo Jo Jo 
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Path integral in G^' can be split into two regions using: 



(Vz) xy = (Vz) xz d 4 z (Vz) zy . (102) 



Therefore the full expression for takes the following form 



f f s OO s OO s OC f 

G (1) = g 2 jd i zjd i z , A(z-z')j dsj rfrj dr' j (Vz) xz (Vz) zz ,(Vz) z , y (103) 



i ft 1 i i-t+t' -i <-s+t+t' 

m 2 (s + r + r') - - / i 2 dr" - - / z 2 dr" - - / z 2 dr" 
4./ 4.L/ 4.L_l^/ 



x exp 

where the intermediate boundary conditions are z = z(t'), and z' = z(t + r'). All 



three path integrals in Eq. 103 can now be integrated to give three free propagators 



between points (x — z — z' — y) 

G (1) = g 2 jd 4 zjd 4 z'A{z-z')S(x,z)S(z,z')S(z',y). (104) 

This is the correct perturbative 1-loop bubble expression, shown in Fig. g^, as 
expected. Next let us consider the tadpole diagram. 



Fig. 22. The simplest bubble diagram 



5.1 Leading order tadpole diagram from the FSR approach 

The tadpole diagram exists only after unquenching. Consider the 1-body prop- 
agator, 



G(y\x) = N\V(f)(detS)S(x,y)e 

In the quenched approximation det(S) is set equal to one. Here let us consider the 
next to leading order contribution to the quenched approximation. We can make 
the following expansion 

detS = exp[tr(logS)] 

[trQogS)] 2 , 



-su 



~ 1 + tr (logS 1 ) + 
where the trace is defined by 

tr(logS) = jd^xixllogSlx) . 
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Note that the trace is infinite and requires regularization, but this does not effect our 
discussion regarding the derivation of the perturbative result. Ordinary propagator 
has a path integral representation, 



S 



T>z exp 



f 



-sm 



dTz 2 (T)—g / dr 4>(z(t)) 



The logarithm of the propagator can also be expressed in the form of a path integral 

o 1 



lOgS: 



T>z exp 



drz\r)-g \ dr </>(z(t)) 



Therefore the diagram with 1-loop connected to the propagator can be written as 

/ OO / OQ ^ / / 

G{x,y) = J dsj j (Vz) xy j (Vzi) exp 



-m 2 s — m 2 si 



f z 2 dr - [ ' z 2 dr + g 2 [ I 
Jo Jo Jo Jo 



drdri A(z(r) - Zi(n)) 



The next to leading order contribution to the 1-loop-l-particle connected propaga- 
tor is given by 



G(x,y) 



/OO /OO j f / (S f S\ 

J dsj Jlj(Vz) xy j(Vzi)j drj d Tl A 



fJ 

x exp 



(z{r) - Zl ( Tl )) 



-m 2 s — m 2 si 



z 2 dr 



z 2 dr 



(105) 



where the loop trajectory z\ is a circular trajectory and therefore has no fixed initial 
or final coordinates. Using the identity Eq. 100 and splitting the path integral as 
before one obtains 



G(x,y) = g* 



/ OO f OO / OC I OO 

I ds dr dri exp 

Jo ) o s l + T i ) o Jo 



— m 2 (s + si + t + ti) 



x I {Vz) xz(r) d A z{r) {Vz) z(T)y | d 4 
xA(z(t) - z;(r ; ))exp 



dr- 



S+T 



dr- 



z? dr — 



zo 

Sl+Ti 



where the boundary conditions are 



z(0) = x, z(t) — z, z(s +r) = y, 
zi(0) = zi(si + T t ) = zo, Zi(n)=zi. 



Let us note that 



exp[— m 2 (s + t)] 



S + T 



dm 2 exp[— m 2 (s + r)]. 



''dr 



(106) 
(107) 



(108) 
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With this replacement all path integrals can be evaluated easily, as done earlier in 
the 1-body self energy calculation, and the 1-loop Greens function to 0(g 2 ) reduces 
to 



G(x,y) = 



'jd 4 zj 



d 4 zi S(x, z)S(z, y)A(z - zi) 



f dm 2 \ 

) m 2 ) 



d 4 z S m (z Ql zi)S m (zi,z ). 



Zq and m 2 integrals can be performed by noting that: 

( d 

d 4 z S m (z ,z)S m (z,z ) = -——?S m (z,z) 



dm 2 



lim S m (x,y) 



0. 



— ►oc 

Therefore the 1-loop diagram to 0(g 2 ) is found as 



G(x,y) 



'jd 4 zjd 4 



zi S(x, z)S(z, y)A(z - zi)S(zi, zj) 



(109) 
(110) 

(111) 



which is the expected result from the perturbation theory, as shown in Fig. |25| . 
This concludes the discussion of perturbative results within the FSR method. It 



Fig. 23. The leading order connected 1-loop diagram 

should be clear that one may extend this discussion to higher order diagrams and 
obtain the correct perturbation theory results. It should be noted that unquenching 
to all orders in the FSR approach is numerically not feasable. This is related to 
the fact that the FSR approach relies on the discretization of trajectories. Every 
loop involved in the calculation represents a new discretized trajectory. Therefore 
inclusion of all loops is practically not possible. However one maybe be able to 
extract information about the effect of unquenching by making an expansion in the 
number of loops, that is by introducing loops order by order. More work needs to 
be done on this topic. 

6 Conclusions 

In these lectures the FSR representation has been introduced with various appli- 
cations to scalar held theories. It has been shown that the FSR is an efficient and 
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rigorous method for doing nonperturbative calculations in field theory. The FSR 
approach uses a covariant path integral representation for the trajectories of parti- 
cles. Reduction of field theoretical path integrals to path integrals involving particle 
trajectories reduces the dimensionality of the problem and the associated compu- 
tational cost. The FSR uses a space-time continuum. There are no boundaries in 
space-time and rotational symmetry is respected. 

Applications of the FSR approach to 1 and 2-body problems in particular shows 
that uncontrolled approximations in field theory may lead to significant deviations 
from the correct result. Results presented here indicate that the ladder approxima- 
tion for the 2-body bound state problem, and the rainbow approximation for the 
1-body problem are both poor approximations. In both cases the crossed diagrams 
(such as crossed ladders) play an essential role. 
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